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We study the two-nucleon force at next-to-next-to-leading order in a chiral effective field theory 
with explicit A degrees of freedom. Fixing the appearing low-energy constants from a next-to- 
leading order calculation of pion-nucleon threshold parameters, we find an improved convergence 
of most peripheral nucleon-nucleon phases compared to the theory with pions and nucleons only. 
In the delta-full theory, the next-to-leading order corrections are dominant in most partial waves 
■ considered. 

(N 

L_i 1 PACS numbers: 13.75.Cs,21.30.-x 



(N 



X 



I. INTRODUCTION 



> : 

J"""-. , The forces between nucleons based on chiral effective field theory have been studied in great detail over the last 
OO • decade, for reviews see e.g. [1, 2]. Most calculations are based on the chiral effective Lagrangian formulated in 
terms of the asymptotically observed ground state fields, the pions and nucleons chirally coupled to external sources. 
The excitation of baryon and meson resonances is encoded in the low-energy constants (LECs) of the pion-nucleon 
interaction beyond leading order. Such a framework provides an accurate representation of the nucleon-nucleon (NN) 
phase shifts if extended to sufficient high order, state-of-the-art calculations have been carried out to next-to-next- 
*"^% ■ to-next-to-leading order in the chiral expansion. This is very similar to the description of pion-nucleon scattering 
in chiral perturbation theory well below the A-excitation energy, where a precise description of the pertinent phase 
shifts has been obtained in a complete one-loop (fourth order) calculation, see e.g. Ref. [3]. Still, it can be argued 
that the explicit inclusion of the delta, which is the most important resonance in nuclear physics, allows one to resum 
a certain class of important contributions and thus leads to an improved convergence as compared to the delta-less 
theory, provided a proper power counting scheme such as the small scale expansion (SSE) [4] is employed. The SSE 
, is a phenomenological extension of chiral perturbation theory in which the delta-nucleon mass splitting is counted 
as an additional small parameter. This improved convergence has been explicitely demonstrated for pion-nucleon 
scattering where the description of the phase shifts at third order in the SSE comes out superior (inferior) to the 
third (fourth) order chiral expansion in the pure pion-nucleon theory [5]. Note, however, that the theory with explicit 
deltas has more LECs at a given order due to the richer operator structure when spin-3/2 fields are present. Clearly, 
for a description of pion-nucleon scattering in the A-region or pion production in NN collisions the inclusion of the 
spin-3/2 fields as active degrees of freedom is mandatory, see e.g. [6] for a review. 

In this work, we want to analyze the two-nucleon forces by systematically including the A-resonance beyond leading 
order in the small scale expansion, extending earlier work presented in Refs. [7, 8] (the precise relation to these papers 
will be given below). More precisely, we construct and analyze the two-pion exchange NN potential at next-to-next- 
to-leading order (NNLO). For that, we must determine the sum of the subleading ttNA LECs 63,6s together with 
the ttN LECs Cj (the values of these differ from the ones obtained in chiral perturbation theory). This is achieved 
by working out the threshold coefficients of pion-nucleon scattering from tree graphs at second order which is of 
sufficient accuracy for the inclusion of the corresponding operators in the NNLO potential. From that, we can deduce 
the momentum and coordinate space representations of various isoscalar and isovector potentials. Having done that, 
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we calculate the peripheral phases in NN scattering based on perturbation theory. It is well established that these 
peripheral waves do not require a non-perturbative resummation and therefore let one most directly analyze the effects 
of various contributions to the effective two-nucleon potential. 

The manuscript is organized as follows. In Sect. II we work out the contributions of the two-pion exchange (TPE) 
graphs including the A- resonance at NNLO. We then determine the corresponding dimension- two LECs from a fit to 
the ttN threshold parameters in Sect. III. The resulting coordinate space representation of the TPE potential and 
the peripheral NN waves are shown and discussed in Sect. IV. We end with a summary and outlook. The appendix 
contains explicit analytical formulae for the ttN threshold coefficients. 



II. A-CONTRIBUTIONS TO THE TWO-NUCLEON FORCE UP TO NNLO 



In this chapter, we construct the two-pion exchange potential (TPEP) with one or two intermediate A-state(s) at 
NNLO. We employ here standard Weinberg power counting for the two-nucleon effective potential, that is the leading 
TPEP starts at next-to-leading order (NLO) (y = 2) and the first corrections to it appear at NNLO (y = 3). The 
leading order (y = 0) potential is given by the static one-pion exchange. Since we are only considering peripheral 
waves with the angular momentum I > 2, no four-nucleon contact interactions contribute to the accuracy we are 
working at. 



A. Effective Lagrangian and power counting 

The calculations performed in the following are based on the effective chiral Lagrangian of pions, nuclcons and deltas. 
We employ here the heavy baryon formulation and display only the terms of relevance for our study: 

C& + £( 2 ) + . . . 

N 4^ N + fAi 1 lT+ (f N + h.c.) , 

NA^ N N+(TA^ NA N + h.c) , (2.1) 
iv ■ D + qau ■ S , 

ci (X+) + c 2 (v ■ u) 2 + c 3 u ■ u + c 4 [S* M , S^u^Uv , 
- (iv ■ D 11 - (m A - mjv) S ij +...) , 

P^ a i(b 3 + b 8 )w i af3 v^ , (2.2) 

where N denotes the large component of the nuclcon field, T is an abbreviation for the large component of the delta 
field, T = T* with i an isospin and /i a Lorentz index. We use standard notation: U(x) = u 2 (x) collects the pion 
fields, — i(v) ' d^u — ud^v)), x+ — u^x u ^ + U X* U includes the explicit chiral symmetry breaking, (. . .) denotes a trace 
in flavor space and D^D 1 ^) is the chiral covariant derivative for the nucleon (delta) fields. Furthermore, is the 
standard projector on the 3/2-components, = g^v— v^v v — 4S , /J S' J/ /(1— d), with the four-velocity, the covariant 
spin vector and d the number of space-time dimensions. We also have w l a = (t 4 u q )/2 and w l afi — {r l [d a , up])/2. In 
the following, we also use the notation A = toa — toat for the TV A mass splitting (which can not be confused with the 
same symbol denoting the delta field). For further notation and discussion, we refer to Ref. [5]. The pertinent LECs 
are at leading order the nucleon axial- vector coupling gA and the 7riVA axial coupling Ha- At NLO, we have the four 
dimension- two ttN LECs a (i = 1,2,3,4) and the combination of 7riVA LECs b 3 + i> 8 . Strictly speaking, all couplings 
and masses appearing in the effective Lagrangian should be taken at their chiral limit values, but to the accuracy we 
are working, we can use their pertinent physical values. 
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FIG. 1: Leading contributions to the 27r-exchange NN potential with single and double A-excitations. Dashed, solid and 
double-solid lines represent pions, nucleons and A-isobars, respectively. Solid dots denote the leading order (LO) vertices. 



The various terms are ordered according to the so-called small scale expansion, in which the small expansion parameter 
Q includes external momenta, pion masses and the nucleon-delta mass splitting, 

Q e {p/A x , M*/A. x , (m A - m N )/A x } , (2.3) 

with A x ~ 1 GeV the scale of chiral symmetry breaking. The various contributions are ordered according to the chiral 
power v > 0. Note also that in the Weinberg power counting used here, 1/ra corrections to vertices and propagators 
are suppressed by an extra power of Q. We remark that while the SSE is a phenomenologically viable extension of 
chiral perturbation theory, it can not be straightforwardly used to study the chiral limit of QCD. For that, one has 
to enforce decoupling. We refrain from a more detailed discussion of this issue here since it does not play a role for 
the results presented in this paper. 



B. Two-pion exchange potential with intermediate deltas 



The NN potential in the center-of-mass system (CMS) can be conveniently expressed in the form: 

V = Vc + Ti • T 2 W C + [Vs + Ti • T 2 W S ] CTi • (72 + [V T + T 1 ■ T 2 W T ] Si ■ q (7 2 • q (2.4) 

+ [V L s + ti • r 2 W LS ] i(a! + a 2 ) -(qxk) + [V aL + t±-t 2 W aL \ ai-{qx k) a 2 ■ (q x k) , 

where the superscripts C, S, T, LS and aL of the scalar functions Vc, ■ ■ ., W„l refer to the central, spin-spin, 
tensor, spin-orbit and quadratic spin-orbit components, respectively. Further, p and p' are the initial and final CMS 
momenta, <7i (r^) refers to the spin (isospin) matrices of the nucleon i and q = p' — p, k = 2 {p' + p)- The leading 
contributions to the 2NF due to intermediate A excitations arise at NLO, v = 2, from diagrams shown in Fig. 1. In 
the context of chiral EFT, they were first discussed by Ordonez et al. [7] using old-fashioned time-ordered perturbation 
theory. These contributions were then re-considered by Kaiser et al. [8] using the Feynman graph technique. In that 
work, compact analytical expressions for the corresponding non-polynomial pieces of the TPEP were presented. For 
the sake of completeness, we list below the results from [8] which are generalized to spectral-function regularization 
with an arbitrary cutoff A (for a precise definition of spectral function regularization, sec [9]): 



• A-excitation in the triangle graphs: 

Wc = ~ 2Tfl% { (6S _ W2)LA((Z) + 12A ^ DA (l)} ■ ( 2 - 5 ) 

• Single A-excitation in the box graphs: 

Vc = S§A {2M " + q2)2AA{qh 

Wc = -^| i {(12A 2 -20M2- llq*)L K {q) +6E 2 Z? A (g)}, 
Vr - -^Vs = -^,{-2L\q) + ^-A^)D\q)Y 
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• Double A-excitation in the box graphs: 

Vc = -^^{-^ 2 LHq) + mH<i) + (z+8A 2 )DHq)}y 

Wc = -i^%{( 12S -^) iA (9) + 3S[^ A ( g ) + (8A 2 -S)^ A (g)]|, 
Vt = -^ = -^%{ 6 ^) + ( 12A2 -- 2 )^ A («)}' 

W T = -^ s = - Y ^^[2LH q ) + (^+^)DH q )). (2.7) 



The quantities S, L A , A A , D A and H A in the above expressions are defined as follows: 
£ = 2M 2 + q 2 - 2A 2 , 

L\q) = 0(1 2M„) !± In ± + 2Ag " S , -=V?T4Mf, S = V^ 7 ^ 
W V J 2q 4M 2 (A 2 + g 2 ) V * ' V 
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The subleading contributions to the 2NF due to intermediate A-excitations arise at NNLO, v — 3, from the diagrams 
shown in Fig. 2. The contribution from the first graph was already discussed in Ref. [7] based on time-ordered pertur- 
bation theory. Notice that the loop integral was treated numerically, and no contributions involving the subleading 
7T-/VA vertices were considered in that work. 

The simplest way to evaluate the subleading A-contributions to the NN potential is using the Feynman graph tech- 
nique. Since the corresponding diagrams do not involve reducible topologies, one obtains the same result for the 
corresponding static TPEP using Feynman diagrams, time-ordered perturbation theory or the method of unitary 
transformations. We find the following expressions for the subleading A-contributions: 



• A-excitation in the triangle graphs: 

Vc = - I ^{6s(4c 1 M 2 -2c 2 A 2 -c 3 (2A 2 + S))^( g ) 

+ ( - 24 Cl M 2 + c 2 {lo 2 - 6E) + 6 c 3 (2A 2 + £)) L A (q) J , 

Wc = {h i + os^ A { 12A2 s ° Hq) + ( -" 2 + 6S)£A(g) } ' 

Wt = ~Wa = {(c 2 - 4A a )D A («) - 2L a (<d) ■ (2.9) 

• Single A-excitation in the box graphs: 

Wc = - \ b i 9 ^ A A DHq) + - 12(A 2 + £))**(<,)} , 
Vt - -±V 8 = - ih \ b ^% hAA {(- 2 - 4A 2 )^) - 2L A (q)} • (2.10) 
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FIG. 2: Subleading contributions to the 27r-exchange NN potential with single and double A-excitations. Filled squares denote 
subleading vertices. For remaining notation see Fig. 1. 



• Double A-excitation in the box graphs: 

Vc = 6W C = ~ 4{b3 ^2 ) J? A { 3 ( 8A2 - S ) + 3SgA (g) + (-^ 2 + 12S)£ A (g)| , 

Vt = ~V S = 6W T = -^W S = ih A {(-c 2 + 12A 2 )^( g ) + 6L A ( g )| . (2.11) 

It is instructive to verify the consistency between the results obtained in chiral EFT with and without explicit 
A's. Clearly, both formulations differ from each other only by the different counting of the AiV mass splitting, 
A <~ M v <C A x versus M x <C A <~ A x . Expanding the various terms in Eqs. (2.5)-(2.11) in powers of 1/A and 
counting A ~ A x should, therefore, yield either terms polynomial in momenta (i.e. contact interactions) or non- 
polynomial contributions absorbable into a redefinition of the LECs in chiral EFT without explicit A's (in harmony 
with the decoupling theorem). Expanding Eqs. (2.5)-(2.11) in powers of 1/A and using the relations 

D Hl) = -^L\ q )-^{AMl + q *)L\ q )+polyn+0(^) , 

H K (q) = L K {q) - L K {q) + polyn + O (A~ 4 ) , (2.12) 

where polyn refers to terms polynomial in q 2 , it is easy to verify that the only remaining non-polynomial contributions 
up to NNLO are the ones given in the first and the last lines in Eq. (2.6). They can be exactly reproduced in EFT 
without explicit A's by an appropriate shift in the LECs C3 and C4 [8]. The contributions to these LECs due to the 
A, C3 = — 2c4 = — 4ft^/(9A), were analyzed already in [10]. 



FIG. 3: Tree graphs for pion-nucleon scattering at NLO. The filled circles/squares denote leading/subleading vertices. Crossed 
graphs are not shown. Note that due to the Weinberg counting of the nucleon mass, some graphs with fixed dimension-two 
LECs do not appear at this order. Also, diagrams that simply lead to a nucleon mass shift are not shown. 



III. DETERMINATION OF THE LECS FROM PION-NUCLEON SCATTERING 



Next, we must fix parameters. At leading order, we use qa = 1-27 for the nucleon axial-vector coupling. The 
corresponding iVA axial-coupling h a is less well known, we therefore use here two extreme values, namely Jia = 1-05 
as in Ref. [5] and Ha — 3(?a/(2v / 2) = 1.34 from SU(4) (or large N c ). To proceed, we must determine the various 
dimension-two LECs. We determine these LECs from a fit at order Q 2 to the S- and P-wave wN threshold parameters, 
which is consistent at the order we are working. Note that because of our treatment of the nucleon mass, there are no 
1/m-corrcctions to the wN amplitude at this order. The pertinent Feynman diagrams describing elastic irN scattering 
at this order are shown in Fig. 3. The explicit analytical expressions for the threshold parameters calculated at 0(Q 2 ) 
are collected in App. A. These depend on the ttN LECs Ci (i = 1, . . . , 4) and the combination of irNA LECs b 3 + b s . 
To determine these couplings, we fit to the EM98 phase shifts, corresponding to Fit 2 in Ref. [5]. The resulting values 
for the LECs are collected Tab. I and the corresponding threshold parameters are shown in Tab. II. Various remarks 
are in order. First, we note that the inclusion of the A leads to much smaller values for c 2 and c 3 and a sizeably 
reduced C4. This is consistent with the findings in [10], where it was shown that about half of the strength of C4 is 
due to the p-meson exchange in the t-channel whereas C2 and C3 are largely given by the delta. The two fits with 
the delta, denoted fit 1 and fit 2 in Tab. I differ by the value of Iia used as input. The resulting values of the LECs 
are quite different, but (as we will show later) the TPEP does not change much. The combination 63 + 6g comes out 
of natural size in both fits. Also notice that both fits lead to very similar results for irN threshold coefficients (the 
differences are not visible in Table II). We also stress that there are significant deviations between different PWAs 
leading to further uncertainties in the LECs. Furthermore, additional significant changes in the LECs appear if one 
goes to 0(Q 3 ). For these reasons, we do not assign any uncertainties to these LECs in our study which can only be 
addressed with sufficient precision in a one-loop calculation within the SSE. What is already obvious from the values 
for the dimension- two LECs is that the TPEP contribution will be greatly enhanced/reduced at NLO/NNLO leading 
to a better convergence than for the delta-less theory. This is also reflected in the results displayed in Tab. II. At 
second order in the chiral expansion one is not able to fit some of the threshold parameters, most notably the P-wave 
scattering volumes. This changes drastically when the delta is included - this is consistent with the statement made 
earlier that the inclusion of the A resums some important contributions to the nN scattering amplitude (but not all). 



IV. RESULTS FOR THE POTENTIAL AND THE PERIPHERAL WAVES 



Having determined the LECs we are now in the position to present the results for the potential and peripheral partial 
waves. To be precise, we first consider the various components of the NN potential in coordinate space and then 
discuss the resulting peripheral waves. 

The coordinate space representations of the various components of the TPEP up to NNLO are defined according to 



V(r) = V c (r)+W c (r) (n ■ t 2 ) + 



V s (r) + W s (r) (n • r 2 ) 



(<7l • <T 2 ) 
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LECs 


Q 2 , no A 


Q 2 , fit 1 


Q 2 , fit 2 


Q A no A [11], fit 2 


Cl 


-0.57 


-0.57 


-0.57 


-1.42 ±0.03 


C2 


2.84 


-0.25 


0.83 


3.13 ±0.04 


C3 


-3.87 


-0.79 


-1.87 


-5.58 ± 0.01 


C4 


2.89 


1.33 


1.87 


3.50 ±0.01 


/lA 




1.34* 


1.05* 




63 + &8 




1.40 


2.95 





TABLE I: Determinations of the LECs from S- and P-wave threshold parameters in irN scattering based on the Q 2 fits with 
and without explicit A's. LECs used as input are marked by the star. Also shown are the values determined in Ref. [11] from 
fit 2 at Q 3 without explicit A's (the errors given are purely statistical and do not reflect the true uncertainty of the LECs). 
The LECs a and bz + bg are given in GeV -1 . 





Q 2 , no A 


Q 2 fits 1, 2 


Q 3 no A [11], fit 2 


EM98 




0.41 


0.41 


0.49 


0.41 ±0.09 


&0 + + 


-4.46 


-4.46 


-5.23 


-4.46 


a o+ 


7.74 


7.74 


7.72 


7.73 ± 0.06 




3.34 


3.34 


1.62 


1.56 


a -_ 


-0.05 


-1.32 


-1.19 


-1.19 ±0.08 


at 


-2.81 


-5.30 


-5.38 


-5.46 ±0.10 


a i+ 


-6.22 


-8.45 


-8.16 


-8.22 ±0.07 


a+ + 


9.68 


12.92 


13.66 


13.13 ±0.13 



TABLE II: Values of the S- and P-wave threshold parameters in -kN scattering for the various fits as described in the text on 
comparison with the data. Both fit 1 and fit 2 lead to very similar results. Units are appropriate inverse powers of the pion 
mass times 10~ 2 . 
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The functions Vc,s,t(t) can be determined for any given r > via 

1 
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2ir 2 r 



2Mr 



d\x\xe » r p c {n) ■ 



V T {r) - 
V s (r) = - 
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where the spectral functions Pi{p) are obtained from Vi(q) by [12]: 

Pi(fi) = Im [V t (0 + - ip, 



dppe ^'(3 + 3pr + p 2 r 2 )p T (p,), 
dfifie~^ r (p 2 p T (p) - 3p s (p)j , 



(4.14) 



(4.15) 



The coordinate space representations of the isovector parts are given by the above equations replacing Vc,s,t(i°) 
Wc\s,r{r) and Vc,s,T(q) — * Wc,s,T(q)- Notice further that for p > 2M,r: 
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FIG. 4: Isoscalar two-pion exchange potentials in coordinate space. Left (right) panels show the results with (without) explicit 
A's. Dashed and solid lines refer to the NLO and NNLO results, respectively. There are no contributions to Vc (Vt,s) a t 
NLO-A (NNLO-A). 
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In Figs. 4 and 5 the results for the functions Vc,s,T(r) and Wc,s,t{t) entering the TPEP are shown up to NNLO in 
EFT with (EFT-A) and without explicit A's (EFT-^k) and using A = 700 MeV. For the LECs Cj and 63 + 6 8 , we use 
the values obtained in fit 1, see Tab. I. The results at NNLO-4& are based on the values for Ci given in the second 
column of Tab. I (these are not very different from the ones obtained in the literature using other values of these 
LECs). The isoscalar central potential Vc(r) is known to be the strongest part of the chiral TPEP. In EFT-4&, it 
receives the first contribution at NNLO. The SSE leads to a much more natural convergence pattern for Vc(f) with the 
dominant contribution being generated at NLO-A. The correction at NNLO-A is still significant and produces further 
attraction. Notice also that the NNLO-4& result for Vc{r) lies between the NLO-A and NNLO-A ones. Contrary to 
the central part, the isoscalar tensor and spin-spin potentials Vr,s(r) as well as the isovector central potential Wc(r) 
receive the contributions at NLO-4& with no further corrections at NNLO-4&. As shown in Figs. 4 and 5, the leading 
contributions due to the intermediate A-excitation turn out to be significant in all these cases while the subleading 
corrections (i.e. at NNLO-A) are less important. Notice that the resulting potential Wc(r) at NNLO-A is nearly 
twice as strong as at NNLO-^L The strength of the remaining isovector tensor and spin-spin parts of the potential 
Wt,s(t) a t NNLO-^k comes out remarkably close to the one obtained at the same order in the EFT without explicit 
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FIG. 5: Isovector two-pion exchange potentials in coordinate space. Left (right) panels show the results with (without) explicit 
A's. Dashed and solid lines refer to the NLO and NNLO results, respectively. There are no contributions to Wc (Wt,s) a t 
NNLO-^ (NLO-A). 



A's, see Fig. 5. Again, contrary to EFT-^k, the NLO contributions to Wt.sM i n the theory with explicit A's do not 
vanish. Let us now comment on the sensitivity of the TPEP to the value of the LEC Iia used as an input in the 
determination of the LECs Cj and 63 + 63, see the discussion in section III. While the values of the LECs obtained in 
fits 1 and 2 for two different choices of Ka differ significantly from each other, see Tab. I, the resulting S- and P-wave 
threshold parameters nearly coincide in both cases, see Tab. II. As demonstrated in Fig. 6, a similar trend holds 
also for the TPEP. The only visible (but still very small) difference is observed for the isoscalar central part of the 
potential. 

We will now consider D- and higher partial waves up to NNLO in chiral EFT with and without explicit deltas following 
the lines of Refs. [8, 9]. Using Born approximation for the scattering amplitude, the phase shifts and mixing angles 
in the convention of Stapp et al. [15] are determined by the 2N potential V as: 

S ' j = -J^WWW ' ^ = - (i^' " l! SjMj + ls SJ) ' (4 - 1?) 

where s and j refer to the total spin and angular momentum, respectively, and q is the nucleon CMS momentum. 
Clearly, such an approximation violates unitarity. This violation is small provided that the corresponding phase shifts 
are small. Alternatively, one can use the X-matrix approach in order to restore unitarity, see e.g. [16]. Here and in 
what follows, we adopt the same notation for the matrix elements in the \lsj) basis as in Ref. [13]. The expressions 
for the partial wave decomposition can be found in appendix B of this reference. 
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FIG. 6: Sensitivity of the NNLO two-pion exchange potential to the choice of the low-energy constant Ka- Solid (dashed) lines 
show the NNLO potentials with Ka = 3<?a/(2V2) ~ 1-34 {Jia = 1.05) and the LECs Cj chosen according to fits 1 and 2 in 
Tab. I. 



Our results for D- F- and G-waves and for the mixing parameters £2.3,4 are shown in Figs. 7-9. The two sets of the 
LECs obtained in fit 1 and fit 2, see Tab. I, yield nearly the same results for all phase shifts and mixing parameters as 
one might expect from Fig. 6. In all channels, the convergence pattern is strongly improved in the SSE as compared to 
the pure chiral expansion. While the NLO corrections are small and the dominant contributions in the EFT-A arise 
at NNLO, the changes between NLO and NNLO are typically much smaller than the ones between LO and NLO in 
the EFT-A. The only exceptions are observed in the 3 -Di, and 3 G§, where the NLO-A result appears to be close 
to the LO one. In all other partial waves, the results at NLO-A are very similar to the ones at NNLO-A. While most 
peripheral partial waves are reasonably well described at both NNLO-A and NNLO-A using the SFR cutoff A = 700 
MeV, the description of the 3 D 3 and 3 Gs phase shifts is rather poor. It should, however, be understood that since 
the phase shifts in these channels are much smaller in magnitude than in the other D- and G-waves, respectively the 
observed deviations have little effect on NN scattering observables. Notice also that, especially in the partial wave, 
the contributions due to iteration of the potential in the Lippmann-Schwinger equation, which are not considered in 
this work, might be significant. 

Finally, we would like to emphasize that the results shown in Figs. 7-9 correspond to one particular choice of the SFR 
cutoff A = 700 MeV. The theoretical uncertainty associated with the variation of A in the range A = 500 . . . 800 MeV 
is exemplified in Fig. 10 where the 3 F± partial wave is shown at NNLO-A and NNLO-A. We note that the uncertainty 
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FIG. 7: D-wave NN phase shifts and the mixing parameter £2- The dotted curve is the LO prediction (i.e. based on the 
pure OPEP). Long-dashed (short-dashed) and solid (dashed-dotted) lines show the NLO and NNLO results with (without) the 
explicit A-contributions and using the SFR with A = 700 MeV. The filled circles (open triangles) depict the results from the 
Nijmegen multi-energy PWA [17, 18] (Vir ginia Tech single— energy PWA [19]). 



at NNLO is comparable in both approaches. Note further that the results for NNLO-4& differs slightly from the one 
displayed in Ref. [9] due to a different choice of the LECs 03,4. For a detailed discussion on the theoretical uncertainty 
in both peripheral and low partial waves in EFT without explicit A's the reader is referred to [9, 13, 14]. 



V. SUMMARY AND CONCLUSIONS 



In this paper, we have analyzed the two-nucleon two-pion-exchange potential in an effective field theory with explicit 
deltas (based on the so-called small scale expansion) at next-to-next-to-leading order in the chiral expansion. The 



12 




50 100 150 200 50 100 150 200 

Lab. Energy [MeV] Lab. Energy [MeV] 




Lab. Energy [MeV] 

FIG. 8: F-wave NN phase shifts and the mixing parameter £3. The dotted curve is the LO prediction (i.e. based on the 
pure OPEP). Long-dashed (short-dashed) and solid (dashed-dotted) lines show the NLO and NNLO results with (without) the 
explicit A-contributions and using the SFR with A = 700 MeV. The filled circles (open triangles) depict the results from the 
Nijmegen multi-energy PWA [17, 18] (Vir ginia Tech single— energy PWA [19]). 



pertinent results can be summarized as follows: 

i) At NNLO, the potential with single and double delta-excitations is shown in Figs. 1,2. It depends on 7 low- 
energy constants, the leading order axial-vector couplings qa and Iia, the dimension-two irN LECs ci, 02,03,04 
and the combination 63 + from the ttNA Lagrangian. 

ii) We have determined these dimension-two LECs from a fit to the pion-nucleon S-and P-wave threshold parame- 
ters, based on the 0(Q 2 ) representation of the ttN scattering amplitude in the SSE, see Tables I, II. As expected, 
the values for the c, (i — 2, 3, 4) are considerably smaller than in chiral perturbation theory. 

iii) The resulting coordinate space potentials show a more natural convergence behaviour than in the theory without 
deltas. In the delta-full theory, the leading contributions to the TPEP at NLO turn out to be significant while 
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FIG. 9: G-wave NN phase shifts and the mixing parameter £4. The dotted curve is the LO prediction (i.e. based on the 
pure OPEP). Long-dashed (short-dashed) and solid (dashed-dotted) lines show the NLO and NNLO results with (without) 
the explicit A-contributions and using the SFR with A = 700 MeV. The filled circles depict the results from the Nijmegen 
multi-energy PWA [17, 18]. 



the subleading NNLO corrections are smaller, see Figs. 4,5. 

iv) This convergence pattern is also reflected in the peripheral partial waves, that have been calculated in first 
Born approximation, see Figs. 7-9. Furthermore, for most partial waves the description at NLO in the delta-full 
theory is better than in the delta-less theory. We also find a mild improvement in some partial waves at NNLO. 
For the larger D-waves, these results should only be considered indicative due to the accuracy of the Born 
approximation. 

v) The theoretical uncertainty due to the variation of the SFR cutoff A is of the expected size and comparable to 
the case of the delta-less theory at NNLO, see Fig. 10 for one typical case. 
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FIG. 10: Theoretical uncertainty in the 3 Fa partial wave due to the variation of the SFR cutoff as described in the text. Left 
panel: Delta-full theory. The upper (lower) band refers to NNLO (NLO) and the dotted line is the leading order OPEP. Right 
panel: Same for the delta-less theory. 



The work presented here paves the way to a systematic analysis of nuclear forces based on a theory with explicit 
deltas. A lot of work remains to be done before this will be achieved, in particular: 



- The low partial waves and the deuteron require a non-perturbative treatment, i.e. the solution of the regularized 
Lippmann-Schwinger equation. This will be dealt with in a subsequent paper [20]. 

- To achieve a higher precision, we have to go to N 3 LO. To do that, one must first reconsider irN scattering 
within the SSE at one- loop order, extending the work of Ref. [5]. It also requires a careful re-evaluation of the 
corresponding two- and three-pion exchange contributions with intermediate deltas. That these corrections will 
not be negligible follows from the calculations of sub-leading three-pion-exchanges in the delta-less theory, see 
Ref. [21]. 

- In the theory with explicit deltas, three-nucleon force (3NF) arises already at NLO. Its leading term is included 
in all current 3NF models and is usually referred to as the Fujita-Miyazawa force [22]. The first corrections to 
it appear at NNLO, some of these are again proportional to the combination of LECs 63 + 63. The ordering of 
various contributions to the 3NF in the theory with explicit deltas is therefore very different from the one based 
on the pure pion-nucleon theory, as discussed e.g. in [23]. 

- Additional isospin violation effects arise in the 2NF and 3NF due to the mass splittings of the various delta 
charge states. These should be investigated systematically within the consistent EFT extending earlier model- 
dependent work [24]. 
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APPENDIX A: tvN THRESHOLD COEFFICIENTS AT ORDER Q 2 

In this appendix we give explicit expressions for irN threshold coefficients a^ + , and af ± at order Q 2 (calculated 
in the SSE) used in section III to determine the corresponding LECs. These are in standard notation: 

< = air^m + M,) (-^+^ + ^)mM^ 

8irF ^m{m + M ) { 2ci( - 2m ~ M -) M - + c 2( 4m " - 2mM - + M l) + c 3( 4 ™ 2 - 2mM « + M l)) > 
1 



a i- = - 



m fcfi - 4c 4 M.) + 27 ^ a(ro + ^ )(Mx + A) - - 2^ + W) , 
- + 2c 3 M.) + 277r ^ 2(m + ^ )( ^ + A) - - 2/M (6a + 6 8 )M.) , 
m (5a + 2c 4 M T ) 

+ M7rJ g (m + M x )(Mj-A») m (^ (2M ^ + A) + ^ ih + h)M ^ + 2A) ) ' 
°k = 24 7 r^M^ + M,) w(g ^ 4C3M - ) 

" 27^( TO + M,)(M2-A 2 ) m 0^ + 2A > + 2 ^ ^ + + A ^) • ^ 



32wF%mM 7r (m 


+ M W ) 


1 




24irF*M v (m 


+ M W ) 


1 




12TrF%M w (m 


+ M w ) 


1 
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